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Abstract 

We recently presented a series of dark energy theorems that place constraints on the equation 
of state of dark energy (wde), the time- variation of Newton's constant (G), and the violation of 
energy conditions in theories with extra dimensions. In this paper, we explore how current and 
future measurements of iude and G can be used to place tight limits on large classes of these 
theories (including some of the most well-motivated examples) independent of the size of the extra 
dimensions. As an example, we show that models with conformally Ricci-flat metrics obeying 
the null energy condition (a common ansatz for Kaluza-Klein and string constructions) are highly 
constrained by current data and may be ruled out entirely by future dark energy and pulsar 
observations. 
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I. INTRODUCTION 



Beginning with the work of Kaluza and Klein [IH3] and continuing today with string 
theory and M-theory, extra dimensions have been a common feature of unified theories. 
The basic notion is that the observed 3+1-dimensional universe is actually described by a 
general relativistic theory in a space-time with one or more extra compactified dimensions. 
If the compactification scale is much greater than 1 TeV (or << 1CT 16 cm), laboratory 
experiments, even at the Large Hadron Collider, are unable to uncover direct evidence of 
extra dimensions. 

In this paper, though, we show how measurements of the equation of state of dark energy 
{w~de) and the time variation of Newton's constant G/G can be used to test or rule out the 
existence of extra dimensions for large classes of models independent of the compactification 
scale. This surprising power to discriminate among extra-dimensional models derives from 
a set of "dark energy theorems," first described in [5H7]. 

The theorems are based on the observation that the expansion of the usual three large 
dimensions tends to cause extra dimensions to vary with time, which, in turn, causes a change 
of u>de and G in the corresponding 4d effective theory. These changes can be avoided in 
a decelerating universe by introducing conventional interactions strong enough to keep the 
sizes of the extra dimensions fixed. However, the dark energy theorems show that, once 
the universe starts to accelerate, conventional interactions satisfying the classical (strong, 
weak and null) energy conditions no longer suffice no matter the size of the extra dimensions. 
Many well-motivated extra-dimensional models satisfy one or more of the energy conditions. 
For these large classes, the dark energy theorems, combined with the observed acceleration 
rate, can be used to compute the predicted time-variation of w^e and G for given model 
parameters. As illustrated below, current measurements are already strong enough to rule 
out a substantial range of model parameters. The more exciting prospect is anticipated 
improvements in the measurement of wde, as described by the Dark Energy Task Force 
(DETF), and of G/G, as constrained by pulsar timing, that can test or rule out whole 
classes of extra-dimensional models. 

The dark energy theorems were derived for the general case of k extra spatial dimensions, 
but the predictions depend on k. For the purposes of illustration, we focus in this paper 
on the well-motivated class of 9+1-dimensional theories (k = 6) with a conformally flat 
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Ricci (CRF) metric and satisfying the null energy condition (NEC) - theories commonly 
used in string- and M-theoretic phenomenological and cosmological models. In Ref [7], we 
showed that this class of models is inconsistent with the standard ACDM model which has 
wde = — 1. In fact, the dark energy theorems show that, the closer wde is to —1, the more 
rapidly the extra-dimensional volume and, hence, G must vary. Furthermore, if the theory 
contains no mechanism for violating the null energy condition, Wde cannot remain close to 
— 1 for an extended period. This means that improved limits on G, combined with ever- 
tightening bounds on the time-variation of wde from future experiments, can progressively 
constrain or rule out this entire class of extra-dimensional models. 

A corollary of this analysis is that a dark energy mission, even if it fails to find any time- 
variation of wde and is consistent with wde = — 1, can still be highly informative because 
it would eliminate well-motivated extra-dimensional models. A second corollary is that a 
coordinated effort is needed. The ambitious improvements in the measurements of wde 
alone, as projected by the DETF, or of G/G alone, as estimated from planned pulsar timing 
surveys, are not sufficient. Each constrains some range of parameter space but leaves some 
substantial untested range. The two approaches are complementary, though: by pursuing 
both to some degree, entire classes of extra-dimensional models can be tested and ruled out. 

This paper is organized as follows. In Section [TTJ, we introduce the conformally Ricci- 
flat (CRF) class of extra-dimensional models used to exemplify our approach and review 
the constraints imposed by the dark energy theorems on the total equation of state w tot 
and Newton's constant G. Here and throughout this paper, the symbol w tot refers to all 
contributions to the energy density of the universe (matter, radiation, etc.), not just dark 
energy. The current value is w to t = —0.74 based on observations [8]. We use the symbol 
wde to refer to the dark energy component alone. 



In Section III , we show how to translate these constraints into predictions for dark energy 



and pulsar timing experiments. The key results are in Section IV, where we compare the 
predictions with current measurements and near-future experiments. In particular, we show 
how current measurements of dark energy and pulsar timing and anticipated improvements 
can be used to test and perhaps rule out the entire class of models. In Section [Vj we conclude 
with a discussion of the generalization to other classes of extra-dimensional models. More 
details on the constraints on the time-variation of G and the dark energy equation of state 
in extra-dimensional models are given in the Appendices. 
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II. DARK ENERGY THEOREMS AND BASIC EQUATIONS 



In this section, we present a brief review of the "dark energy theorems" first described and 
proven in [SHTj. The theorems impose constraints on extra-dimensional theories in which 
the four- dimensional effective theory undergoes cosmic acceleration (w to t < — 1/3). The 
theorems assume a spacetime metric of the form 

ds 2 = 9MN dX M dX N = e 2n ^g™ w (t, x) dxW + h a p{t, y) dy a dy? (1) 

Here, /i, v ... are indices along the four large spacetime dimensions with coordinates x M , and 
a, f3 ... indices along the k compact dimensions with coordinates y a , and N,M ... assume 
values along both compact and noncompact dimensions. In (fxl) , is a flat Friedmann- 

Robertson- Walker metric, which we take to be 

where m,n range over 1...3. The extra-dimensional metric h a p and warp factor Q can be 
time-dependent. For the purposes of illustration, we restrict ourselves in this paper to 
"conformally Ricci-flat" (CRF) metrics: 

K P it J y)=e- 2 ^hll{t J y) (3) 

where h^J(t, y) has vanishing Ricci scalar, and Q is the same warp factor which appears in 
0. The CRF metric occurs in many string theory models such as warped Calabi-Yau [13] 
and warped conifold jS] constructions (where they are sometimes referred to as conformally 
Calabi-Yau metrics). 

In this work, we further assume that the extra-dimensional matter satisfies the null energy 
condition (NEC): 

T MN n M n N > (4) 

for any null vector n , where Tmn is the stress-energy tensor. For perfect fluids in four 
dimensions, the NEC requires that p + P > 0. The NEC is satisfied by scalar fields 
with canonical kinetic terms (regardless of their potential), de Sitter and anti-de Sitter 
cosmological constants, p-form fields, and positive-tension extended objects. To violate the 
NEC requires exotic ingredients such as scalar fields with higher-derivative kinetic energy 
terms or negative-tension extended objects (such as orientifold planes). Very often violation 
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of the NEC leads to ghosts, instabilities, problems with gravitational thermodynamics, and 
other pathologies. 

We will refer to the spectrum of extra-dimensional models having CRF metrics and 
obeying the NEC as the NEC/ CRF family of models. For concreteness, we will take the 
number of extra dimensions to be k = 6. 

As shown in [5H7j, constraints can be found by considering extra-dimensional theories 
for which the 4d effective theory is described by Einstein gravity and has w tot < —1/3, 
and, then, following the overall expansion and dilation of the extra-dimensional metric in 
the extra-dimensional Einstein equations. The time evolution of the higher-dimensional 
metric can be expressed asthe combination of a trace part £ and a symmetric, traceless 
shear component a aj 3 [5H7]: 

where h a ^o a ^ = and h a p, a a p and £ are functions of t and y. The variable £ is the local 
expansion rate of the extra-dimensional space. To measure the overall dilation of the extra 
dimensions, we define a variable (a by 

Ca=^ je A ^Vhd k y (6) 

where H is the four-dimensional Hubble rate, and A is a constant which may be chosen 
for convenience. Hence (a represents the fractional growth of the extra-dimensional volume 
per Hubble time, computed using an A-dependent measure. For the choice A = 2, the 
volume measure in ^ matches the one which determines the four-dimensional Planck mass 
in warped compactifications. Hence, for this value of A we have 

f - -*C (7) 

where G is the four-dimensional Newton's constant, and we dropped the subscript ( 2 = C 
for economy of notation. When ( is nonzero, the volume of the extra dimensions is changing 
with time, and, hence, the four-dimensional Newton's constant is changing as well. 

The dark energy theorems are derived by assuming that the higher-dimensional matter 
fields satisfy the NEC (5H7J. By dividing the space-space components of the stress-energy 
tensor into two blocks corresponding the non-compact and compact directions, two pressure- 
like parameters can be constructed by taking trace averages over the 3x3 and k x k blocks 
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of the higher-dimensional metric, 

p 3 = l -g mn T mn and p k = \h^T afi , (8) 

where m, n range over the spatial coordinates 1...3 and a, /3 range over the extra-dimensional 
coordinates as in 0. The NEC is violated if either p + p 3 or p + pk is less than zero at any 
space-time point, where p = — T ° is the higher dimensional energy density; or if the volume 
weighted average of either is less than zero; or if either of the "A-weighted" averages (the 
volume- weighted averages of e An (p + p 3 ) or e An (p + pk)) is less than zero for any A. By 
combining the higher-dimensional Einstein equations, expressions can be derived relating 
the A-weighted averages to £, cr Q( g, w tot , fl, k and the 4d effective energy density p^. The 
condition that the A = 2- weighted average of p + Pk be non- negative can be rearranged into 
the constraint: 

> a + at( + a 2 ( 2 (9) 



dN 

where iV = ln(a) and a is the Einstein frame scale factor. The analogous condition for p + p 3 
is 

C 2 < F (10) 

The functions a , «i, «2, and F depend on w tot , the number of extra dimensions k and 
the type of extra-dimensional metric. For our example with CRF metric and k = 6 extra 
dimensions, the strongest observational constraints are found for A = 2, in which case the 
value of £ is related to the variation of Newton's constant in four dimensions by (J7|). Then, 
the constraint equations ^ and (10) required for the extra-dimensional theory to satisfy 



the Einstein equations and the NEC become [5H7]: 



and 



d C > , 2 , 3(w to t ~l) c 9(1 + 3wto t ) 



e < 9(1+Wtot) . (12) 



With (11) and (12) in hand, we can derive some qualitative consequences of the dark energy 
theorems for the NEC/CRF family of models. For example, suppose we wish to find a 
solution for which the four-dimensional Newton's constant does not vary (£ = 0). With this 



choice of £, the constraint equation (11) can only be satisfied if w to t > —1/3. Conversely, 
if the four-dimensional universe is accelerating, then the ( must be non-zero and the 4d 
Newton's constant must vary. 



If Wtot < — 1/3 and the universe is accelerating, ( cannot vanish but one might look for 
cases where w to t and ( are constant. Such solutions require that there exists a value of ( 



such that the right-hand side of ( 11 ) vanishes. This is only possible when u> tot > —5 + 2-^/5 ~ 
—0.53. We will denote this value of w to t by u>trans- Thus, a second corollary of the dark 
energy theorems is that steadily accelerating solutions with constant ( are only possible 
when wtot > w tmns , which in this case corresponds to w to t ~ —0.53. 

The most striking conclusions of all are reached by considering values of w tot below 



Wtrans- F° r "^tot < ^trans, the right-hand side of (11) is positive definite for all values of 



C, so ( must evolve with time. Since ( is increasing with time, after a finite number of 



e-foldings the bound encapsulated in (11) and (12) will be violated. Hence there are no 
steadily accelerating solutions for u> to t < w tmns ~ —0.53: acceleration with w to t < w tr ans is 
necessarily transient. 

III. OBSERVATIONAL CONSTRAINTS AND THE w - w a PLANE 



The two constraint equations Eq. (11) and (12) restrict the time-variation of G and 
w to t for any extra-dimensional model described by a CRF metric and obeying the NEC. In 
this section, we explain how to express these theoretical constraints as limits on the time- 
dependence of wde and, then, how to incorporate observational constraints that can further 
restrict or perhaps rule out the remaining, theoretically allowed possibilities. 

The first step is to parameterize the time-dependence of wbe and (. In general, each can 
have complicated time-dependence. However, in keeping with standard practice, we take 
wde to be a simple function of the scale factor a over a range encompassing the present 
epoch (0 < a < 2 where a = 1 today): 

wde(a) = w + (1 - a)w a + (1 - a) 2 w b . (13) 

The quadratic form is used, rather than the linear DETF parameterization, because the 
analysis would produce artificially strong constraints for a pure linear dependence. 

For (, there is only one free parameter that, without loss of generality, can be chosen to be 
(acc, the value when a = a acc . For a given {u> , w a , Wb, (acc}, the behavior of ( is determined 
by integrating (|9]); see Ref. [6] for details. We define the set of "theoretically allowed" 
values of {w , w a , w b , (acc} to be those that satisfy all of the following conditions: 
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wde > — 1 for all a < 2, so that the dark energy component satisfies the null energy 
condition for all times; 

wde < as a — > 0, so the universe is sure to be matter- and radiation-dominated at 
very early times; 



• the constraint equations (11) and (12) are obeyed, so the theory is compatible with 
the dark energy theorems. 

• Qde, the ratio of the dark energy density to the critical density equals 0.74, consistent 
with current observational constraints for wde near —1. 

We then project this set of points in the four- dimensional parameter space into the two- 
dimensional wo — w a subspace. A point {wo,w a } in this subspace is labeled "compatible 
with NEC and CRF" if there is at least one choice of {wb, Cacc} such that {wq, w a , Wf,, Cacc} 
satisfies all the conditions above. As shorthand, the set of compatible points is labeled C. 
See Figure [TJ 

Note that the NEC/CRF compatible region in Figure [l] includes the "de Sitter" point 
(wo,w a ) = (—1,0), which appears to contradict the claim in Section [II] that a universe with 
constant ti>tot < Wtrans violates the dark energy theorems. There is no real inconsistency, 
though: a true 4d de Sitter universe with w to t = — 1 (or, equivalently, wo = —1, w a = Wb = 0) 
for all time is incompatible with the dark energy theorems. However, it is also possible to 
choose Wb non-zero such that (wo,w a ) is transiently equal to (—1,0) in the present epoch, 
and yet have to tot increase to values greater than w tra ns in the past and future, in keeping 
with the dark energy theorems. 

We next add the observational limits on the time- variation of wde and G. For wde, the 
current limits on wq and w a and the DETF projections for future experiments are used, as 
reviewed in Appendix |Bj The current variation of G can be computed from the value of ( 
by 

£ = (-7.4 x 10- 11 yr- 1 ) x Co, (14) 

^ today 

where (o denotes the present value of (; the Hubble parameter is taken to be H = 72 
km/s/Mpc; and we have used ^ to relate the variation in G to £■ The present-day instan- 
taneous constraint on G/G is 

| = (0 ± 5) x KT 12 yr- 1 (15) 
8 
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FIG. 1: The (wq, u> a )-plane where wde = WQ+w a (l — a)+Wb(l — a) 2 . The plot shows the boundary 
of C, the set of points compatible NEC and CRF for the case of k = 6 extra dimensions. For any 
point {wo,w a } lying within C, one can find at least one choice of {uibXacc} that satisfies the NEC 
and the dark energy theorems and that is matter-dominated in the past. Conversely, for points 
outside C, at least one of these conditions is violated for every choice of {wf,, ( aC c}- 



which results from various experimental studies described in Appendix |Aj Integrating (J7|) 
gives a formula that relates C(N) to the secular variation of G between two times ti and to 

ni+ \ rN(to) 
K l) - exp / C(N) dN (16) 



For the secular variation of G, the current bound is 

= 1-001°;?°, (17) 

where Gbbn and Go are the values of G at big bang nucleosynthesis (BBN) and at the 
present, respectively. Further details about the constraints on secular variation of G may 
be found in Appendix [A] 
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Using these constraints, we can define a x 2 measure of the agreement between observations 
and models compatible with the dark energy theorems. The four parameters wq, w a , u>& and 
( acc suffice to predict the time variation of u>de and G/G. Then, x 2 as a function of these 
four variables is 

2( r v ^ P +i y ( Wa \ 2 ( g/g y ^ wgp-i y ^ 

X (w ,w a ,w 6 ,Ca CC ) = + t — + rr • — + -=fpj J7TT 5 (18) 



<5u> p 7 \ 5u>„ 



V <5[G/G] J + V 5[Gbbjv/G 



where 5wo, Sw a , ... are la uncertainties, and w p is the value of u>de at the "pivot" red shift 
z p = 0.3. We introduce w p because the DETF convention expresses experimental constraints 
in terms of Wq and w p , as described in Appendix [Bj In terms of our model parameters, w p 
is given approximately by 

w p = w + 0.231w a + 0.0533w fe (19) 

In our definition of x 2 , we have assumed that future observational data returns a result that 
best fits (wq, w a ) = (—1, 0) in the WQ-w a plane, and null results for instantaneous and secular 
variation of G, albeit with some uncertainty in these values. We assume the future data 
differs from current data only in its progressively tighter error bars. 



We have found it useful to reduce the full four-variable x function (18) to one which 



depends on only two variables: Wq and w a . We define this two-variable x 2 function as the 



minimum value of the full x 2 function (18) over all (wb,(acc) in the compatibilty region C, 
for the given (wo,w a ). Formally, 

X 2 (w ,w a ) = min x 2 (w , w a , w b , (acc), (20) 

where the minimization is understood to be over all (wb, (acc) such that (wq, w a , Wb, (acc) £ C. 
If we were estimating the values of model parameters, it would be more appropriate to 
marginalize over u>& and £ acc than to seek the minimum x 2 ■ However, we want to determine 
whether any model compatible with the dark energy theorems can also be consistent with 
future high-precision experiments which give results supporting a cosmological constant and 
with no time- variation in G. Since we are not interested in the specific value of any of our 
model parameters, but instead in the quality of the best fit, x 2 minimization is more ap- 
propriate than marginalization. In the following section, we use the methodology described 
here to compute x 2 { w 0i w a) based on current and near-future observational constraints. 
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IV. RESULTS 



In the previous section, the dark energy theorems were shown to forbid equation of state 
parameters (w , w a ) that lie outside region C in Fig. [T] In this section, we consider how 
observational constraints on wbe and G/G can further limit and perhaps rule out region C 
itself and, hence, the entire family of NEC/CRF models. Current data is consistent with 
u>de = — 1 and G/G = 0; as noted in the previous section, we assume for the purposes 
of this study that improved measurements will continue to point to the same conclusions 
about wde and G/G in the future, but with smaller uncertainties. We compute the minimum 



X 2 (wq, w a ) statistic in (18) based on this assumption. A model will be considered "ruled out 
if x 2 ( w o, w a) exceeds 3a (99.7% for two parameters). 

Figure [2] illustrates x 2 ( w o, , w a ) in region C based on current observations. Only a small 
sliver of C is ruled out. Figure [3] shows that improving measurements of dark energy only 
or of G/G only is not powerful enough to rule out the entire NEC/CRF family of models. 
For example, the figure shows that a substantial range of C near the "de Sitter" point 
(w ,w a ) = (—1,0) is still allowed even if the projected sensitivity of the most ambitious 
and optimistic DETF space-based proposal is achieved. Similarly, a tenfold improvement in 
G/G limits from pulsars, with no dark energy information included, leaves unconstrained a 
substantial range of C with wq ^ —0.9. The key point, though, is that the poorly constrained 
regions of C for the two measurements do not overlap, suggesting that a combination of the 
two can be effective in ruling out all of the NEC/CRF family of models. 

For example, constraining NEC/CRF models by improving dark energy constraints only 
and combining them with current G/G constraints can rule out the entire NEC/CRF family 
of models, at the cost of a very ambitious dark energy mission. Figure [4] shows x 2 { w ®i w a) 
assuming the most ambitious and optimistic DETF concept, a space-based Stage IV proposal. 
All of C is ruled out by more than 3a. (Stage III and Stage IV ground-based missions may 
be able to rule out the entire plane if the optimistic DETF projections hold true (at ~ 3. la 
and ~ 3.4a, respectively), but they would fall short for DETF pessimistic projections.) 
Although the figure is not shown here, improving constraints on G/G only and combining 
with the current limits on u>de is insufficient to rule out all of C. In short, improving only 
one of the two measurements is a difficult approach, at best, for ruling out the NEC/CRF 
family of models. 
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FIG. 2: Current constraints on the (wo,w a ) plane for NEC/CRF family of models with k = 6 
extra dimensions. This plot combines present constraints on both dark energy parameters and G, 
derived using the techniques described in the manuscript. The contours are generated by choosing 
the model which best agrees with experimental constraints amongst all models which obey the 
conditions of the dark energy theorems at each value of (wo,w a ). 

A less demanding strategy involves modest improvements to measurements of both u>de 
and G/G, as illustrated in Figure [5] This figure shows the exclusion regions with DETF 
Stage II iode measurements, and only a factor of two improvement in the current value of 
G/G. The entire range of C can be ruled at ^ 3.4<r. 

V. DISCUSSION 

The results of the previous section demonstrate that it is possible to test and possibly rule 
out an entire class of extra-dimensional models in a way that does not depend on the size of 
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FIG. 3: Imposing the constraints from dark energy or G/G independently does not constrain the 
NEC/CRF family of models. Even a very ambitious dark energy measurement, or a significant 
improvement in pulsar constraints, would leave a region of parameter space which is compatible 
with the NEC/CRF assumption. Hence, we cannot rule out this family of models by considering 
these constraints independently. 

the extra dimensions. The approach relies on the fact that the expansion of the universe is 
accelerating today and that acceleration causes the volume of the compactified dimensions 
to change in cases where the theory satisfies a classical energy condition. Instead of probing 
the compactified dimensions directly, the approach is to test for the effect of time-variation 
of the compactified dimensions on G and u>de- A key advantage of this approach is that the 
time variation required by dark energy theorems does not depend on the size of the extra 
dimensions. 

For the purposes of illustration, we have focused on CRF models that obey the NEC 
because they are common to many string- and M-theoretic constructions. We have demon- 
strated that current observations allow a significant range of these models, but that improve- 
ments in the measurement of G/G and wde anticipated over the next few years can rule 
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FIG. 4: Projected constraints after a Stage IV space-based mission (optimistic) for models obeying 
the NEC with CRF-type metrics in k = 6 extra dimensions. With dark energy measurements alone, 
very ambitious experiments are needed to exclude the full NEC-compatible region, even when 
current constraints on G/G are included. The Stage IV space-based mission is the only DETF 
scenario which excludes the entire plane to > 3<r for both pessimistic and optimistic projections. 



them out entirely. To accomplish the task by improved dark energy constraints alone would 
require the DETF's most ambitious Stage IV plan combined with current measurements of 
G/G. 

The most cost-effective approach,, though, is to combine a modest improvement in G/G 
and wde limits, which can impose constraints tighter than those obtained from a very high- 
precision dark energy measurement alone. Fig. [6] summarizes how different combinations 
of improved measurements of G/G from pulsar timing measurements and dark energy ex- 
periments described by the DETF can rule out the entire range of k = 6-dimensional CRF 
models obeying the NEC at the 3a level or higher. 

A similar approach can be applied to test other metrics and/or other energy conditions. 
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FIG. 5: Projected constraints after a DETF Stage II dark energy measurement and a twofold 
improvement in G/G bounds, for models obeying the NEC with CRF-type metrics in k = 6 extra 
dimensions. This family of models is excluded at the same level as with a Stage IV dark energy 
measurement and current G/G bounds, hence modest improvements in constraints on both dark 
energy and G/G can do much better than ambitious improvements in dark energy constraints 
alone. 




For example, we have carried out a similar analysis for the family of models in which the 
extra-dimensional theory satisfies the strong energy condition (SEC) and is described by 
curved metrics of the form: 

ds 2 = e 2n(t,y )g FRW^ %) + ^ y) ^ 

where g™ w is a flat Friedmann-Robertson- Walker metric; and the extra-dimensional metric 



h a /3 and warp factor f2 can be time-dependent. Unlike the CRF family, in (21) we allow 
the extra-dimensional metric h a p to have arbitrary Ricci curvature. This family of models 
satisfies a set of dark energy theorems that is different from the theorems for NEC/CRF 
family. We find that, as in the case of the NEC/CRF metrics, the SEC/Curved family 
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FIG. 6: The improvements in G/G constraints required to rule out NEC-satisfying models with 
CRF-type metrics and k = 6 extra dimensions, as a function of the available dark energy con- 
straints. For each DETF stage, the bars show the improvement in pulsar constraints that would 
lead to exclusion of this family of models at the 3a, Aa, and 5<r levels. By contrast, this family can- 
not be ruled out (even at the la level) from dark energy information alone, even with constraints 
from extremely ambitious measurements. 



of models is not ruled out by current constraints on w DE and G/G, but DETF Stage II 
limits combined with current constraints on G/G is sufficient to rule out the entire family 
of models. 

We expect to be able to extend our approach to yet more combinations of metrics and 
energy conditions. Even for the NEC/CRF family of models, it may be possible to derive 
additional dark energy theorems. The current theorems specify conditions that are necessary 
to have cosmic acceleration and still satisfy the energy conditions, but they are not sufficient. 
The current theorems were selected because they are the simplest to prove, as discussed in 
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Ref. [7]. However, there may be stronger theorems that combine with observational con- 
straints of u>de and G/G to produce much more stringent constraints on extra-dimensional 
models. 



Appendix A: Constraints on time-variation of G 



In Section III we described how the higher- dimensional Einstein equations and energy 



conditions are distilled into the expressions rt9|) and (10) for £ = G/(GH). Since rt9j) and 



(10) are inequalities, they cannot be used to predict a unique value for (, but they can place 
tight constraint the range of values which ( can assume. This in turn places constraints 
on the allowed variation of G. The basic technique is illustrated in Figure [7} The plot 
assumes a cosmological model in which wde = —0.8 and is constant in time. The upper 
and lower black curves in Figure [7] are ±y/~F. If the extra-dimensional model obeys the 
conditions of the theorems, then ( must remain in the region between the two outer curves, 



according to (10). This immediately constrains the present value of £, and hence G/G. For 
the NEC/CRF family of models considered in this paper, F = 9(l + Wtot)/2 today is of order 
unity, and so this initial constraint tells us that G should vary by no more than order unity 
per Hubble time. 

To determine the constraints derived from the differential inequality ([9]), it suffices to 
consider trajectories (± that saturate the inequality, as illustrated in Figure [7j Since this 
differential equation is first-order, the trajectories obtained with different initial values of ( 
never cross each other, so trajectories that satisfy the inequalities are bounded by (±. As 
a practical matter, we impose initial conditions for ( at the beginning of the accelerating 
epoch, when w = —1/3, and integrate a = 2, figuring that our simple parameterization of 
dark energy should not be trusted for a much wider range of a. 

We define as the solution to ^ which saturates the inequality and uses the smallest 
allowed initial value of ( when ti>tot = —1/3. By an "allowed" initial value, we mean one for 
which the trajectory resulting from integrating ^ never leaves the envelope ±y/~F before 
the end of the integration interval. In some cases, such as the one illustrated in Figure [7], 
the initial value for (_ is at the lower edge of the ±vF envelope, at ( = —yF. In other 
cases, a larger initial value may be necessary to avoid leaving the ±\/F envelope before the 



end of the integration interval. The function Q+ is defined in a parallel fashion: it is the 
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w = -0.8 w a = w b =0 




FIG. 7: The constraints on £ imposed by the dark energy theorem inequalities (|9j) and (10) for 
a flat matter-quintessence model with fidark = 0.74 and u>de = —0.8 being constant, for which 
t^tot passes through —1/3 at z = 0.78. The curves are computed for NEC-satisfying models with 
CRF-type metrics and k = 6 extra dimensions. The upper and lower curves are the limits ±\/~F 
obtained from (10). The curve (_ is the solution obtained by saturating ^ and using the most 



negative initial value for ( allowed by (10). The curve C+ is obtained in the same way, but uses 



the largest possible initial value for £ which does not cross over \/F by a = 2. All £(a) trajectories 
which are compatible with NEC must lie within the central band. 



trajectory defined by (J9J) which uses the largest possible initial value of £ and does not leave 
the ±a/F envelope before the end of the integration interval. 

Using the two functions (± we can put much tighter constraints on the allowed behavior 
of (. The solutions (± divide the allowed region between ±a/F into three subregions: 

The —\/F < £ < C_ region is forbidden by the energy condition. We can show this by 
contradiction. Suppose there is a curve ((a) which reaches a point (G,ai) located in this 
region but satisfies the conditions of the theorems. Then at the beginning of acceleration, 



18 



C > —yF by (10). By continuity, the curve ((a) must have crossed the curve (- at some 
a* < ai, and in particular gone from above the curve to below. This is a contraction, for at 
the crossing, the slope of £ is less than that of (, and since ( = ( at the crossing point, ([9J 
could not have been satisfied. A similar argument applies to ( + < ( < +\/F. Therefore any 



function ( which satisfies both inequalities (|9j) and (10) is confined to the region < ( < £ + . 
Using the relationships and (16), this constraints on £ translate into constraints on the 
allowed values of the ratio between G at the beginning of the accelerating epoch and today, 
and the present value of G/G. 

For example, in the case shown in Figure [7j the present values of £_ and ( + today are 
approximately 5.87 x 10~ 3 and 2.26 x 10 _1 , corresponding to 



- 1.68 x 10~ n yr" 1 < 



G 
G 



< -4.33 x 10~ 13 yr" 1 



(Al) 



today 



where the upper limit comes from £_ and the lower limit from ( + . This entire range is 



essentially consistent with the current experimental limits (15) at roughly 3a. There is a 
gap between the upper limit and zero, which means that G = constant is not consistent 
with the NEC/CRF family of models. This gap could be explored at 3a with a significant 



improvement of the instantaneous G/G constraints (15) by roughly a factor of ~ 35. The 
limits from secular change in G over the accelerating epoch are 



0.678 < < 0.929 

Go 



(A2) 



where the lower limit corresponds to ( + and the upper one to 



1. Observational constraints on the instantaneous variation of G/G 

There is a long history of searches for variation of Newton's constant G (for a recent 
review, see [llj). In recent years interest in measuring the variation of G (as well as other 
fundamental constants) has been driven by ideas in higher-dimensional unification and mod- 
ified theories of gravity which may explain cosmic acceleration. 

To be useful in our analysis, it is important that we use constraints on the time- variation of 
G that do not make specific assumptions about the functional form of the time-dependence. 
For example, some constraints on the instantaneous variation of G/G are obtained by as- 
suming G(t) ~ t 13 . For some special functional forms, such as this, one can obtain very 
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strong observational limits on G/G by applying constraints on the value of G in distant 
past, {e.g., from big bang nucleosynthesis). We cannot use such constraints, because there 
is no guarantee that the variation of G predicted by the dark energy theorems will follow 
any particular simple form. Similarly, we cannot use constraints based on stellar evolution 
that average over one or more Gyrs to place a limit on the current instantaneous variation 
of G since the current value of G/G may be very different from the average over the last 
Gyr. 

We also must avoid constraints which depend on a particular modification of gravity. By 
assuming a particular gravitational theory, it is sometimes possible to obtain tight constraints 
on theory parameters by other measurements, and then use these parameters to derive a 
constraint on the current value of G/G. One common example is to assume Brans-Dicke 
theory: constraining the Brans-Dicke parameter c^bd by, for example, measuring the Shapiro 
delay can be used to provide an indirect constraint on the current variation of G. We must 
discard such constraints because the assumption of a particular gravitational theory is always 
in the background. 

The ideal constraints on the current variation of G are those which use data over very 
short timescales - such as a few decades - that are extremely small in comparison with 
cosmological timescales, and make minimal assumptions about the gravitational theory being 
tested. No test can ever be completely model-independent, but assuming Einstein gravity 
is a conservative assumption which leaves us at little risk for being misled. 

Solar system constraints are ideal for bounding the variation of G, due to the short 
measurement timescales and nearly Newtonian regime. An early constraint on the variation 
of G from analyzing ranging data to the Viking Mars lander [12] gave a constraint of G/G = 
(2 ±4) x 10~ 12 yr -1 The dominant source of uncertainty was the modeling of asteroid effects 
on the orbit of Mars. Using the same data set, but a different asteroid model, others have 
obtained slightly weaker limits of | G/G \ < 3X10" 11 yr" 1 [T3] and \G/G\ < lxlO" 11 yr" 1 [H], 
so it is possible that the initial uncertainties may have been underestimated. Measurements 
using planets in the inner solar system are less sensitive to uncertainties due to the asteroid 
belt. An early constraint of \G/G\ < 4 x 10~ 10 yr -1 [15] was obtained from radar ranging of 
Mercury. This has been sucessively refined over the intervening decades [ToT - fTS] . Currently, 
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the constraint 

G 

— = (0 ± 2) x 1(T 12 yr" 1 (A3) 
G 

follows from an analysis of a combination of ranging data from Mariner 10, Mercury and 
Venus P2]. 

Pulsar timing measurements have also obtained very stringent bounds on the instanta- 
neous variation of G. In contrast with solar-system tests, post-Newtonian effects in pulsar 
systems, such as the influence of gravitational binding energy and damping due to the emis- 
sion of gravitational radiation, cannot be neglected. There are also theoretical uncertainties 
involving the composition of the bodies in the pulsar system. Usually in these cases one 
must take a phenomenological approach to obtain reasonably model-independent bounds. 
One approach relates the anomalous pulse period derivative 8P to the variation of G by 
SP/P = —2G/G. When applied to the Hulse- Taylor pulsar PSR 1913+16, after successive 
refinements [2"0Tl2"2"] this gives the bound 

% = (4 ± 5) x 1(T 12 yr- 1 (A4) 

A similar analysis of the pulsar PSR B1855+09 gives the slightly weaker bound G/G = 
(—9 ± 18) x 1CT 12 yr -1 . This bound is somewhat more conservative than the one for PSR 
1913+16, since the companion is not a neutron star, hence there are fewer composition- 
dependent uncertainties. 

For our work we adopt the canonical present-day constraint of 

G 

— = (0 ± 5) x 1(T 12 yr" 1 (this work) (A5) 
G 



This uncertainty envelops both the recent constraints (A4) from PSR 1913+16 and the inner 



solar-system tests (A3). This is a conservative constraint, for its uncertainty agrees with the 
pulsar constraint, which is the larger of the two. Furthermore, the mean value is taken to be 
zero, in agreement with the solar-system tests and with the simplest theoretical hypothesis 
that G/G = today. 



2. Observational constraints on the secular variation of G 



Another set of constraints on variation of G arise from the integrated variation of G over 
cosmic history. Many of the same caveats for constraints on the instantaneous variation of 
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G apply here as well: use constraints which do not assume a specific functional form for 
the variation of G with time and which do not make strong model-dependent theoretical 
assumptions. These constraints bound the ratio between G , the value of G today, and the 
value of G at an earlier epoch of cosmic history. 

The most precise constraints of the variation of G with cosmic time have come from big 
bang nucleosynthesis (BBN). Changing the value of G changes the Hubble expansion rate 
during BBN, which affects freeze-out temperatures and the rate of various nuclear reactions. 
This leads to different predictions for primordial abundances of D, 3 He, 4 He, and 7 Li. These 
predictions also depend on the baryon-to-photon ratio r] and the number of relativistic 
species during BBN, which is usually parameterized by the number of light neutrino species 
N u . Allowing (G,T),N U ) to vary led to the early bound of 0.7 > Gbbn/Gq > 1.4 [25] . Later, 
using more recent evidence that N u = 3 and an independent determination of r] from cosmic 
microwave background measurements, this bound was refined to [23] 

^ = (A6) 
In this work we take the essentially identical bound 

GbBN i nn+ 0.20 



G 



1.001 ale (this work) (A7) 



which differs from (A6) by recalibrating the mean value of the ratio to be unity, in accord 
with the simplest theoretical hypothesis that G has been constant since BBN. 

There is an additional assumption which is specific to our analysis. As we describe in 



more detail in Section III, the dark energy theorems bound the variation of G from the 
beginning of the accelerating epoch, where G = G acc to the present day, where G = G$. 
There are no precise measurements of G at the transition from deceleration to acceleration, 
so we will use the ratio of Gbbn/Gq as a stand-in for the ratio Gbbn/Gq. 



Appendix B: Constraints on wde 

The discovery that the present universe is accelerating [25j I2S] has triggered theoretical 
efforts to understand this fact, and experimental searches to better characterize the prop- 
erties of the dark energy that is presumably responsible (for reviews, see [2TH3D] ) . One 
simple possibility, consistent with all current data, is that the dark energy is a cosmological 
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constant. It is also possible that dark energy is dynamical, so that its equation-of-state 
parameter wde varies with time. Present and future experimental data can constrain this 
possibility. A significant problem with applying these constraints is that they usually depend 
on a specific choice for parameterizing the variation of wbe with time. 

An issue in combining observational constraints and constraints from the dark energy 
theorems is how to parameterize the time-variation of u>de- Some studies of supernova data 
[321 [33] choose a parameterization, first suggested in [31J, of the form 

wbe = w + w'z (Bl) 

where wo is the dark energy equation-of-state parameter today, and w' = dw/dz\ z =o- With 
this parameterization, the constraints Wq = — 1.3llo28 anc ^ w ' = l-^-o.to were obtained 
|33j . with the errors in each parameter strongly correlated. If u>de is assumed to be con- 
stant, then the constraints u>de = —1-021 019 has been obtained [33]. Another choice of 
parameterization, proposed in [3H [35] , is 

wbe = w + w a (l - a) = w + w a — — (B2) 

X -(- z 

which assumes dark energy evolution linear in the scale factor a. This has been used in 
some recent supernova analyses [361 EZ] , along with data from cosmic microwave background 
(CMB) and baryon acoustic oscillations (BAO). In Ref. [37] this led to correlated errors on 
Wq and w a of roughly ±0.2 and ±0.7, respectively. In the same reference, by a flat universe 
with -u; DE constant gives constraints of u>de = —1-001 ± 0.071(stat) ± 0.081(sys). 

Analysis of the seven-year Wilkinson Microwave Anisotropy Probe (WMAP) data has also 
given constraints on dynamical dark energy [8]. If the equation-of-state u>de is constant, 
then CMB data alone constrains vj-de = —1.10 ± 0.14 at lcr. Allowing the wde to vary 
linearly with a gives the constraints = —0.93 ± 0.13 and w a = — 0.411q^ at la. Similar 
results are obtained from the WMAP five-year analysis [3H] • 

In this work, we draw from these various analyses, and take today's value il>o of the 
equation-of-state to be 

w = -1.00 ±0.13 (this work) (B3) 

The uncertainty is the same as the WMAP seven-year analysis, which compares favorably 
with the the uncertainties in the supernovae analyses. While the WMAP analysis prefers a 
central value of which is roughly la from Wq — l, we have chosen to fix the central value 
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at —1 for the purposes of our analysis, based on the notion that a cosmological constant 
is consistent with all data sets and is the simplest assumption. For the current constraints 
only, we have chosen to place no limits on w a , since the variety of paremeterizations used in 
the literature make it difficult to compare this parameter between analyses. 

For projections of future dark energy parameter uncertainties, the Dark Energy Task 
Force (DETF) report (39] is used. This report introduces the simple linear parameterization 



in (B2). However, projected uncertainties are given in terms of two different parameters 
(w p ,w a ). The parameter w p is the value of »de at the "pivot redshift" z p m 0.3. The pivot 
redshift is the redshift at which a specific experiment can determine the value of u>de(%) 
with minimum uncertainty. 

The DETF report defines the potential progression in our observational knowledge in 
terms of four stages: 

• Stage I. Current experiments. 

• Stage II. Ongoing experiments related to dark energy. 

• Stage III. Near-term, currently proposed projects. Ref. [39] considers ground-based 
surveys for BAO, cluster lensing, supernovae, and weak lensing. 

• Stage IV. Ambitious long-term projects, such as the Large Survey Telescope (LST), 
Square Kilometer Array (SKA), or Joint Dark Energy Mission (JDEM) 

Projections for both pessimistic and optimistic error ellipses in the (w p , w a ) plane are given 
for Stages II-IV. Furthermore, Stage IV projections are given for both ground-only and 
space-based Stage IV campaigns. 

A summary of the various data scenarios we consider is given in Table [Tj 
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TABLE I: Data scenarios for dark energy measurements. These are taken from the DETF report. 
We show the scenario name (used in our computer code) , the uncertainty Aw p in wde a the pivot 
redshift, and the uncertainty Aw a in w a at the pivot redshift. 
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0.36 
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0.016 
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0.15 


DETF Stage IV space, pessimistic 
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DETF Stage IV space, optimistic 



[1] Nordstrom, Gunnar. Physikalische Zeitschrift 15 504506 (1914). 

[2] Kaluza, Theodor. Sitzungsber. Preuss. Akad. Wiss. Berlin. (Math. Phys.) 966972 (1921) 
[3] Klein, Oskar. Zeitschrift fr Physik a Hadrons and Nuclei 37 895906. (1926). 
[4] S. W. Hawking and G. F. R. Ellis. The large-scale structure of space-time. Cambridge Univer- 
sity Press, 1973. 



[5] D. H. Wesley, arXiv:0802.2106 [hep-th]. 



[6] D. H. Wesley, JCAP 0901, 041 (2009) |arXiv:0802.3214 [hep-th]]. 



[7] P. J. Steinhardt and D. Wesley, Phys. Rev. D 79, 104026 (2009) |arXiv:08 11.1614 [hep-th]]. 



[8] E. Komatsu et al, |arXiv:1001.4538| [astro-ph.CO], 

[9] G. W. Gibbons, in F. del Aguila, J. A. de Azcarraga, L.E. Ibanez (eds), Super symmetry, 
supergravity, and related topics. World Scientific, Singapore, 1985. 
[10] J. M. Maldacena and C. Nunez, Int. J. Mod. Phys. A 16, 822 (2001) |arXiv:hep-th/0007018| . 
[11] J. P. Uzan, Rev. Mod. Phys. 75, 403 (2003) [arXiv:hep-ph/0205340l . 
[12] R. W. Hellings et. al. Phys. Rev. Lett. 51, 1609 (1983) 



25 



[13] R. W. Hellings et. al. Philo. Trans. R. Soc. London A 310, 227 (1983) 

[14] J. F. Chandler, R. D. Reasenberg, I. I. Shaprio et. al. Bull. Amer. Astro. Soc. 25, 1255 (1993) 
[15] I. I. Shapiro, W. B. Smith, M. B. Ash, R. P. Ingalls and G. H. Pettengill, Phys. Rev. Lett. 
26, 27 (1971). 

[16] R. D. Reasenberg and I. I. Shapiro in Atomic Masses and Fundamental Constants, Vol. 5, 

edited by J. H. Sanders adn A. H. Wapstra (Plenum, New York 1976) p 643. 
[17] R. D. Reasenberg and I. I. Shapiro in On the Measurements of Cosmological Variations of the 

Gravitational Constant,edited by L. Halperin (University of Florida, Gainsville FL), p. 71 
[18] I. I. Shapiro in General Relativity and Gravitation, edited by N. Ashby, D. F. Bartlett, and 

W. Wyss (Cambridge University Press, Cambridge, England). 
[19] J. D. Anderson et. al in Proceedings of the 6th Marcel Grossman Meeting on General Relativity, 

Kyoto, June 1991. Edited by H. Sato and T. Nakamura (World Scientific, Singapore), p. 353 
[20] T. Damour, G. W. Gibbons and J. H. Taylor. Phys. Rev. Lett. 61, 1151 (1988). 
[21] T. Damour and J. H. Taylor, Astrophys. J. 366, 501 (1991). 
[22] V. M. Kaspi, J. H. Taylor and M. F. Ryba, Astrophys. J. 428, 713 (1994). 
[23] F. S. Accetta, L. M. Krauss and P. Romanelli, Phys. Lett. B 248, 146 (1990). 
[24] C. J. Copi, A. N. Davis and L. M. Krauss, Phys. Rev. Lett. 92, 171301 (2004) | arXiv:astro- 

|ph/031l"334] . 

[25] A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998) 

|arXiv:astro-ph/980'5201~] . 
[26] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 565 

(1999) [aTxiv" : astro- ph/9812 133]. 
[27] S. Weinberg, Rev. Mod. Phys. 61 (1989) 1. 

[28] T. Padmanabhan, Phys. Rept. 380, 235 (2003) [arXiv:hep-th /02 1 2290] . 
[29] P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003) [arXiv:astro-ph/0207347] . 
[30] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006) |arXiv:hep- 
|th/0603057| . 

[31] A. R. Cooray and D. Huterer, Astrophys. J. 513, L95 (1999) [arXiv:astro-ph/9901097l . 
[32] E. Di Pietro and J. F. Claeskens, Mon. Not. Roy. Astron. Soc. 341, 1299 (2003) [ arXiv:astro- 
ph/0207332]. 

[33] A. G. Riess et al. [Supernova Search Team Collaboration], Astrophys. J. 607, 665 (2004) 



26 



[arXiv: astro-ph/04025 1 2] . 
[34] M. Chevallier and D. Polarski, Int. J. Mod. Phys. D 10, 213 (2001) | arXiv:gr-qc/0009008| . 
[35] E. V. Linder, Phys. Rev. Lett. 90, 091301 (2003) [arXiv:astro-ph/0208512| . 
[36] A. G. Riess et al, Astrophys. J. 659, 98 (2007) | arXiv:astro-ph/0611572| . 
[37] M. Kowalski et al, |arXiv:0804.4142| [astro-ph], 

[38] E. Komatsu et al [WMAP Collaboration], |arXiv:0803.0547| [astro-ph], 
[39] A. Albrecht et al, |arXiv:astro-ph/060959~l~j 

[40] D. Huterer and M. S. Turner, Phys. Rev. D 64, 123527 (2001) [arXiv:astro-ph/0012510| . 
[41] I. Maor, R. Brustein, J. McMahon and P. J. Steinhardt, Phys. Rev. D 65, 123003 (2002) 

[arXiv:astro-ph/0112526] . 
[42] W. Hu and B. Jain, Phys. Rev. D 70, 043009 (2004) | jarXiv:astro-ph/0312395l . 
[43] S. B. Giddings, S. Kachru and J. Polchinski, Phys. Rev. D 66, 106006 (2002) |arXiv:hep- 

th/0105097 . 

[44] I. R. Klebanov and M. J. Strassler, JHEP 0008, 052 (2000) [arXiv:hep-th/0007191| . 



27 



